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Abstract. We present a simple exactly solvable extension of of the Jaynes-Cummings 
model by adding dissipation. This is done such that the total number of excitations is 
conserved. The Liouville operator in the resulting master equation can be reduced to 
blocks of 4 x 4 matrices. 



PACS numbers: 42.50.Pq 
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1. Introduction 

In the field of quantum optics it is typical to deal with open systems in which some 
degrees of freedom are treated as an environment. The simplest example is a two level 
atom in thermal equilibrium with a continuum of modes This accounts for an 
effective spontaneous decay from an excited state to the ground state. To deal with 
this kind of problems, one usually works under a Markovian assumption (environment 
without memory) and as a result one can work with a master equation in the Lindblad 
form [2 J , which also follows from the pioneering work of Kossakowski and collaborators 

[SHU- 

Solutions to this kind of equation, except for a two level atom, involve complicated 
expressions due to the extended space in which one is working. The master equation 
can be written as the action of a super-operator or Liouville operator £ on a density 
matrix p. If one has annxn density matrix, the super-operator will act on a n x n 
dimensional space. 

In this work we present a new form for a Jaynes Cummings model with dissipation. 
The aim here is simply exact solvability with Kraus operators not considered previously. 
We show how to construct a master equation whose Liouvillian can be diagonalized in 
blocks of 4 x 4 matrices and thus produce analytical solutions. Yet we will see that in 
the same way more complicated problems, such as spins larger than 1/2 can be brought 
to finite matrix block form and therefore at least allow numerical treatment to arbitrary 
accuracy. 

Indeed in two recent papers we displayed models solvable in closed form for two 
different atoms in a cavity and for Dirac-Moshinsky oscillators coupled to an isospin field 
[5j E] . These models were based on the fact, that Hamiltonians could be constructed, 
such that the total number of excitations was an additional conserved quantum number. 
In the basis where this quantity was diagonal the Hamiltonian broke into blocks no larger 
than 4x4. 

Here we make an analogue construction for the master equation in the Lindblad 
form for the Jaynes Cummings model in rotating wave approximation using Kraus 
operators that preserve the total number of excitations. This Hamiltonian reduces to 
2x2 blocks and the corresponding super-operator to 4 x 4 blocks again thus guaranteeing 
solvability. 

In section[2]we present a review of the Jaynes- Cummings model, which will allow us 
to fix the notation to be used in section[3]where we introduce the explicit master equation 
that describes a system with dissipation, but at constant number of total excitations. 
This allows a spectral decomposition of the Liouville operator in terms of 4 x 4 matrices. 
In section [4] we present explicit solutions to the eigenvalue problem when the atom is in 
resonance with the cavity and discuss the dynamics of two particular initial states. 
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Consider the Jaynes-Cummings Hamiltonian in the rotating wave approximation 

H = 5cr z + g (ao~+ + aV_) . (1) 

One can identify an additional conserved quantity, which can be interpreted as the 
number of excitations 

J = aW £ (** + !)• (2) 

The basis in which / is diagonal is given by the states 

|n-l,l), \n,0), (3) 

where |n — 1, 1) represents a state with n — 1 photons in the field mode and an excited 
atom. |n, 0) accounts for n photons and the atom in the ground state. The action of / 
on each of the previous states is given by 

I\n- j,j) = n\n- j = 0, 1. (4) 

In this basis H is block-diagonal, with the blocks 



( ' *f , (5) 



with the corresponding eigenenergies 



±E n = ±y/8*+g*n. (6) 

The eigenstates -the dressed states- of H are given by 

= cos# n |n - 1, 1) + sin0 n |n, 0), 
\4>~) = - sin^n - 1, 1) + cos^ n |n, 0) (7) 

with 



6 n = arctan J — ^ . (8) 




3. Master equation 

The Heisenberg equation for a density matrix in a closed system is given by 

p=-i[H,p\ (9) 

For an open system description one can consider a master equation in the Lindblad form 
which describes a non- unitary Markovian evolution [2]. It can be written in terms of an 
arbitrary set of Kraus operators Oj as 

p = Cp 

= {0]0 3 p + pOp, - 20 3 pO\) . (10) 
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Where we have introduced the Liouville super-operator C that acts on a density matrix 
p. £ is a linear operator and the generator of a completely positive dynamical semigroup 
[3]. Known solvable models [7J El |9] for this type of equation are, for instance, taking 
one operator O = a for a cavity with losses, or O = a z for an atom with spontaneous 
decay. 

In this work we explore the possibility of having other types of operators which 
allow closed solution of the master equation. The condition we impose is that such 
operators commute with J, this means that this dissipation will conserve the number 
of excitations. There are many other operators that one could consider for instance / 
itself, or a^a and a z . Here we arbitrarily chose the pair 0\ = aa + and O2 = aV_. 
With this we can construct the following master equation which describes dissipative 
dynamics 

p = £p= -i [H, p] 

— (aa + a) o~-p + paa + a) a- — 2a)a-pao^j 

— (a) a-acr + p + pa)a-ao + — 2aa + pa)a^j . (11) 

As mentioned above, the operators we have chosen preserve the number of 
excitations /. Then we find it convenient to work in the basis in which I is diagonal to 
to represent any density matrix as 



P= E E p£ m \n-3,j){™-k,k\ = J2pn,m- (12) 

n,m=0 j,k=0 n,m 

Here we have partitioned the density matrix into the 2x2 matrices p n .m with matrix 
elements p J n k m - Note that for n = there is a single state as n — j > 0, which in this 
case implies j = 0. Each p n , m has a definite number of left and right excitations, which 
can be summarized by its commutation relation with / as 

[I,Pn,m] = {n- m)p njm . (13) 

By construction the Liouvillian preserves the number of excitations and maps any 
Pn^m into another p' n m , i.e. it does not couple blocks. Actually one can write an effective 
Liouville equation for each subspace spanned by a pair of excitations n and m as 

Pn,m £"n,mPn,mi (l^) 

where each C n ^ m is a superoperator acting on a 4— dimensional space. The full Liouvillian 
is simply the sum of all these terms and can be expressed as 



If we express each block of the density matrix as a column vector 

\ 



Pn,r, 



rn,m 



\ rn,m 



p l,0 



0,0 
rn,m 



pL,l 
rn,m 

pi.o 

p 0,l 
rn,m 

0,0 
rn.m 



(16) 
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r 



-y/rnn'jo 
igy/m 

-igy/n 



%g\jm 
-2i5 1 



... v 'run j —igy/n 

2iS — \\Jmn^ igy/m 

—igy/n ig\/m —y/mnji 



(17) 



where we have defined 



7i - To- 



lls) 



For m = n each block reduces to the Liouville super-operator of a two-level atom in 
thermal equilibrium with a reservoir and driven by a classical field. 

The problem is now reduced to solving the eigenvalue problem for each non- 
hermitian 4x4 matrix of the Liouvillian. In this case, right and left eigenvectors 
are not the same and the eigenvalue equations read 

r — \U) tfj) 

*-"n,mfn,m n,mrn,m 

xO) r = A (i) (19) 
The biorthonormality condition is 



(20) 



The evolution of an initial density matrix p(0) under a time independent Liouvillian 



reads 



pit) = e a p(0). 



Using the spectral decomposition and biorthonormality this can be rewritten as 

oo 4 



n,m=0 j=l 



with the coefficients 



,0") 



=Tr(p (i) on). 

,m \rn,m rU J 



(21) 



(22) 



(23) 



Now that we have found the formal solution to the problem, we can study the dynamics 
through the evaluation of relevant quantities. Here we shall concentrate on the atomic 
system and consider its reduced density matrix, which can be evaluated as 



n 

to obtain the 2x2 matrix 



(24) 



\ Qoi Qoo J 



( p 1 ' 1 



E 

n 



1,0 \ 
Pn+l,n 



„0,1 n,0 
\ Hn,n+1 rn,n 



(25) 
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One can now evaluate the population inversion of the atomic system 

W{t) = g n (t) - g 00 (t) = 2g n (t) - 1 (26) 

as well as its purity 

P(t) = gUt) + g 2 00 (t) + 2\g 01 (t)\ 2 . (27) 

It is equally interesting to study the behaviour of the cavity Then the Wigner 
or Husimi functions should be evaluated to get some insight of the action of the 
Kraus operators proposed by analyzing the corresponding phase space along the lines 
of ref. [TH]. We leave this for future investigation. 



4. Explicit solutions for zero detuning 



We now concentrate on the case when the atom is in resonance with the mode, i. e. 
zero detuning 5 = 0. Here one is able to find simple explicit solutions for eigenvalues 



and eigenvectors of the Liouville operator of the master Equation (11) by diagonalizing 



the blocks (17). We shall rescale to unity coupling between atom and cavity, that is 
g = 1. The four eigenvalues for each block are given by 

3a 



A (2) 



A (3) 



i ran 



-7 - 



'mn 



4 

Run 



-7 + 



-7 - 



7 2 — ( \/m + 



imn 
1~6~ 



7^ 



m + 



-7 2 - [y/m - y/n 



A (4) 



'mn 



-7 + 



imn 



7 2 



(2* 



Note that for m = n the fourth eigenvalue is zero, which corresponds to the stationary 
state of the dynamics. 

It is convenient to introduce 

= \f™ ~ , X (j) 
n,m 2 ' m,m 

To simplify the notation in the following equations we shall write Aj 
With this notation, the full set of left and right eigenvectors written in matrix form: 



(29) 



A CO and / ■ = 1® 



Pn,m 



~2 

Pn,m 



/ • (2+70^1 )(sM+Jn) 
' 1 l 2 (4+hl) 



wyo+m-yi — 2li 
h^+hl) 
_2±7oh_ 
(4+«27) 



—I 



rrryo+Tvyi — 2l\ 

■ (2+7i<i)(y^7i+Vw) 

m,7o+n7i-2?2 \ 
( v / m+ v / n)(4+i 2 7) 
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Pn,m 



Pn,m 



Pn,m 



-2 

Pn,m 



P 3 

rn,m 



Pn,m 



A3 — A4 A3— A4 



V 



a 3 

A3— A4 
a 3 

a 3 — A4 



sfm-yfn 
A3— A4 



\fm-\Jn 
A3— A4 



A3— A4 



A3 
A 3 — A4 



h-h 

h 
h-h 

h_ 

h-h 



h-h 



—l- 

—i 



( 

«2-d 

\ h—h h—h 
I ^ (2-7iA 4 )(v^-y^) 



h-h / 
/m+y/n \ 
h-h 



h-h 
71 A 

A 3 (4-A 4 7) 



n7o+m7i+2A4 \ 
A 3 (4-A 4 7) 



m7o+«7i+2A4 
A 3 (4-A 4 7) 

2-71 A3 
(4-A37) 



, • m7o+rt7i+2A3 
\ ^(v / ^-v / S)(4-A 3 7) 



■ (2-7 A 4 )( v ^n- v / n) 
1 A 3 (4-A 4 7) 

■ n7o+r»7i+ 2 A3 \ 
i ( v ^- v ^)(4-A 3 7) 



2 -70 A3 
4-A37- 



(30) 



(31) 



For m = n the proper limit has to be taken in the fourth eigenvector and that is 

/ 4+ra77i „• 2^(70-71) \ 



"4 _ 
Pn,n 



V 



4+ra77i 
8+n7 2 

; 2Vn(7o-7i) 
' 8+nr 2 



8+n7 2 

4+71770 
8+n7 2 



(32) 



Having found explicit solutions for the eigenvalue problem, we are going to use 
them to investigate the dynamics of two simple initial conditions. 



4-1. A dressed state as an initial condition 

In this subsection we investigate the behaviour of an eigenvector of H, which is a 
stationary state if the dissipation is turned off. In this case with zero detuning (5 = 0) 
such state has the form 

|0+) = -L(|n-l,l) + KO>). (33) 

The eigenstates of H have a definite number of excitations n, and thus the initial state 
can be expressed as 



p(O)= Pn , n (O) = |0+)(0+| = 2 
As different blocks do not couple we obtain 



1/11 



1 1 



p(t) = Pn , n (t) = ETr{p£|0+)(0+|}e^pg) 
i=i 



(34) 



(35) 
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Figure 1. Population inversion W (top) and purity P (bottom) as a function of time 
for an initial dressed state. We show 4 curves, all with the same number of excitations 
n = 2 and 70 = 71 for the gray curve, 70 = 0.08, 71 = for the solid black curve, 
70 = 1.2, 71 = for the dashed curve and 70 = 0, 71 = 1.2 for the dotted curve. 



In this case we find that the relevant matrix elements are given by 



i,i /.x _ 2 - 7i AS 70 -71 A l$ n e 

PnnK 1 ) A . (3) _ + ,(1) 7 (2) 



2 \(J) t 



'i An,n J ln,n ln,n j=l o ~r £ln,n J 
1.2 ■ V '"V /" 1 ' 



— / \ A (3) / 

^(70 — TlJ „ ne A n,n l 

1 



e = I T 1 1,1 • (37) 



,. ^(70-71) ^ ^ f _ iy , 

6 7 (i) /(2 ) 2^. ,(,•)- v X J W 

Once we have the density matrix, we can compute the reduced density matrix for the 
atomic system from Equation (25). The off-diagonal terms vanish and we are left with 
the diagonal reduced density matrix for the atom 

To visualize the dynamics we choose to evaluate the population inversion W(t) of the 
atomic state, that can be obtained from Equation (26) and its purity P(t) from Equation 
(27). This can be trivially achieved, noting that for the reduced density matrix (37) the 
matrix element Qu = p^ n gives us all the information about the atomic subsystem. In 
figure [T] we have plotted these quantities, for different choices of 71 and 70. 

First we will discuss the case 70 = 7i- This is plotted in gray and the result are 
the constants values W = and P = 1/2 in figure [TJ Of course this case includes 
the absence of dissipation, where the dressed state is an eigenvector and therefore a 
stationary state. For non zero but equal dissipation constants, we see that the state 
remains stationary. 

For small dissipation constants where we set 70 = 0.08 and 71 = (black curve) we 
see small oscillations that eventually decay to the initial values of P and W. 
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Increasing in an asymmetric way the dissipation constants (70 = 0, 71 = 1.2) the 
oscillations are strongly damped and leading to a greater final population of the excited 
state (dotted line). We find similar results when (70 = 1.2, 71 = 0), but here the ground 
state population dominates for large times (dashed line). For the last two conditions, 
the behaviour of purity is equivalent and leads to a final steady state which is not a 
complete mixture as was the case for the initial state, i.e purity of the atom increases 
due to the external coupling. 



4-2. A superposition with two different number of excitations 

As another example we now consider the time evolution of an initial product state of 
the field and atom subsystems, that is 

|^/o) = cosa|ri, 0) + sina|n, 1), (38) 

The first and second terms in the last expression correspond, respectively, to states of 
n and n + 1 excitations. Taking the outer product of this vector we obtain the initial 
density matrix which results in a mixture of terms with these two excitations. The 
resulting matrix can be expressed as the sum of the four terms 

p(t) = Pn,n{t) + Pn,n+l\t) 1 Pn+l,n\t) + Pn+l,n+l 

(t), (39) 

and can be represented as a 4 x 4 matrix. Every single term can be computed using the 



spectral decomposition as in Equation (22). Taking partial trace over the cavity's degree 
of freedom to obtain the reduced density matrix of the atomic system, as in Equation 
d25k, yields 



Pn+1 n+1 Pn+1, n(t) 



g(t) 



V P°ni+l(t) Pn'%(t)+Pnll,n+l{t) ) 

The 2x2 density matrix is determined by two of its elements as 



Qn(t) = cos 2 a I 



+ 



®(2+©7i: 



(40) 



A (2) t 



+ (4 + (ffl - ln% (4 + l^j) (® - l n % , 



4 + (n + 1)771 4+i,n+i(2 + C+i,n+iTo)e A »+ 1 .»+ 



+ sin a\ 



(i) 



+ (n + 1)7 2 (4 + l 



(i) 

n+l,n+l 



mi 2) 



+l,n+l 



(1) 

n+1, n+1 j 



i 2 i,n + i(2 + /^i, n+ i7o)e A "+ 1 ."+ 



7(2) 



(2) 



+ 



(A _L /( 2 ) ^\(J^ —1 



Qoi(t) 



%+l,n+ll)v'n+\.ii + \ 

cos a sin a 



(i) 

n+1; 



(n 7 + 7i - 2£ +1 )e A ».'« + x 4 (7*7 + 71 - 2l^ n+l )e x ^ 



i 1 ' ,t 



(2) 



(4+£ + iT)a 



(2) _ ,(1) ^ 
n,n+l 'n,n+ly 



(4 + 41+17) (J 



(2) _ ,(1) * 
n,n+l 'n.n+1/ 



+ 



(n 7 + 7o + 2Ai 4 i +1 )e A --+^ (n 7 + 7o + 2Ai 3 i +1 )e A -™+ 



(3) 



(4 - A 



(-1) 



v( 4 ) 



n,n+l 7 )(^n,n !. 



^n!n+l) (4 — A^ +l7 )(A 



(3) 



n,n+l 



A 



(3) > 
n,n+l 1 



(41) 
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Figure 2. Population inversion W (top) and purity P (bottom) as a function of time 
for an initial dressed state. We show 4 curves, all with the same number of excitations 
n = 2 and 70 = 71 = for the gray curve, 70 = 0.08, 71 = for the solid black curve, 
70 = 1.2, 71 = for the dashed curve and 70 = 0, 71 = 1.2 for the dotted curve. 



In the same way as in subsection 4.1[ from the reduced density matrix of the atom, 
one can evaluate the inversion of population W{t) (26) and the purity P{t) (27). In 
figure [2] we have plotted the corresponding values for different choices of 71 and 70 . 

For the gray curves (71 = 70 = 0) we can see that the initial state is not a stationary 
state of the system without dissipation. 

Turning on one of the dissipation constants (70 = 0.08, 71 = 0), the black curve 
shows how the oscillations are damped to end up with a state with almost equal 
population between excited and ground state and close to the total mixture. 

Increasing in an asymmetric way the dissipation constants results in a similar 
behaviour as in the case where we used an initial dressed state. The oscillations are 
strongly damped to end up with a greater occupation of the excited state for for the 
dotted curve (70 = 0, 71 = 1.2), and the opposite for the dashed curve (70 = 1.2, 
7i = 0). 



5. Conclusions 



We have displayed a new open system that can be solved analytically. From the point 
of view of solvable systems this is always an interesting step, both formally and also as 
a test ground for approximate or numerical solutions for problems that are not solvable 
analytically. Solvability here is based on conservation of the excitation number both in 
the unitary and the non-unitary evolution. 

As far as the behaviour of the population inversion is concerned, we note that the 
two terms treat the cavity mode and spin asymmetrically if the two gammas are not 
equal. Thus the observed result is not entirely surprising, but it may well deserve a more 
detailed study because purity will reach its lower limit only with symmetric dissipation 

While this model depends on the exact solvability of fourth order polynomials, since 
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Galois we know that such a general solution does not exist for higher degrees. We could 
try to cook up models that correspond to special solvable cases, but we would rather 
consider the following. If e.g. we have a three or more level atom (or a spin larger than 
1/2) or even say two two-level systems, according to the discussion in the introduction 
the dimension of the matrix representation of the super operator is the square of the 
dimension of the Hilbert space. Thus for spin 1 we would have a 9 x 9 matrix and for 
two spin 1/2 particles a 16 x 16 matrix form of the super-operator. Obviously these 
are not in general diagonalizable in closed form. Yet numerical diagonalization can be 
considered exact in the sense that arbitrary exactitude can be reached. From that point 
on eigenvalues and the dual sets of eigenfunctions can be used to obtain the results 
required for a study of decoherence. 

Other more complicated states, such as coherent states or thermal states can be 
used as initial states, and will be considered in future work. The question how to 
construct a specific experiment that corresponds to these conditions, unfortunately, is 
still an open problem, but we are striving to find a solution. 
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